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To measure the energy spread of a energy electron beam up to 5 GeV in a single shoot 
measurement the possible applications of an existing magnet as spectrometer has been 
discussed. For distinguishing the energy resolution and an optimal experimental setup 
numerical simulation, based on the measured magnetic field, have been performed. 



1 Problem 

Electron pulses, which are aimed on to be studied, 
may have a comparative narrow energy speared (such 
as ~ 5%) pQ [SJ. On the other hand the electrons 
energy may vary between different electron bunches. 
Since the interest lies in knowing the energy distribu- 
tion of a single bunch and not in the average distri- 
bution of the bunches it is required to measure the 
energy spread in a single shoot procedure. In addi- 
tion the small spread in energy requires a high energy 
resolution. 

One possibility to determine the energy distribu- 
tion in a single shoot measurement is given by using a 
dipole magnet as spectrometer. The discussed spec- 
trometer can cover a large energy range (~ lOOMeV 
to 5000Mey) simultaneously. 

Since the desired magnet has been used for a dif- 
ferent purpose (switching magnet) it is important to 
know its magnetic properties. Therefore a map of the 
magnetic field has been generated for different coil 
current. In addition the ramping characteristic of the 
magnet has been tested. 

To investigate the resolution the electron's trajec- 
tory has been studied by numerical simulation. These 
simulation allows to take the real magnetic field into 
account. 



2 Coordinate system 

Depending on the problem it can be more convenient 
to deal with problems either in a lab's rest frame or 
in the beam's frame. In the lab's rest frame two dif- 



ferent coordinate systems are used. The first is the 
cartesian system (x, y, z)), where the z axis is per- 
pendicular to the pole surface, the x axis is along the 
connection between the source and the magnet's cen- 
tre and the y is perpendicular to the x and z axis (Fig. 
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Figure 1: Sketch of the orientation for the different 
coordinate systems. The magnet's pole are pictured 
blue. The beam's trajectory is red and the three co- 
ordinate systems are black. 



The coordinate system (X, Y, Z) is parametrised 
so that the Z axis is along the propagation of the 
beam's centre, the Y axis is perpendicular to Z and 
z. And the X axis is perpendicular to Y and Z (Fig. 
QJ. Due to the radial symmetry of the magnetic field 
it is some time more convenient to use a cylindrical 
coordinate system (r, </?, Z). 
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3 Discussion on the Magnet 

3.1 Magnet's physical data 

The studied magnet was previously used as switching 
magnet. It was fabricated in 1989 by Danfysik (serial 
number 89374). It generates a held of 1.15 T at 87 A 
The physical radius of the magnet is ~ 256 mm the 
magnetic radius is ~ (167.0 ± 0.6) mm (see Fig. [5]). 




Figure 2: Sketch of the magnet which will be used as 
the spectrometer. The hrst picture shows the ground 
plan and the second the sheer plane of the magnet. 
The yoke is blue, the magnet poles are red and the 
coil is green pictured. In the second drawing the coil 
is omitted. 




Figure 3: Photo of the magnet with the old vacuum 
chamber. To measure the B z field the hall probe was 
inserted in the different tubes. 

Since the largest part of the field is contributed 
by the component along z a good knowledge of 
B z is crucial for determining the beam's deflec- 
tion angle. Therefore a profile of the vertical 
field was measured as a function of the radius for 
-30°, -15°, 0°, 15°, 30° and 45° deflection angle. 
These angles were given by the alignment of the 
previous vacuum tube (Fig. [3]), which had ports at 
0°, ±15°, ±30° and ±45°. To insert the hall probe 
into the 

different vacuum tubes the probe was fixed on 
a pipe of stainless steel which fitted in the vacuum 
tubes. The position of the hall probe was controlled 
by measuring the probe's position relative to the vac- 
uum tube's entrance (± 0.5mm). The twisting angle 
f3 between the hall probe and the field was maximally 
0.5°. If the B field is assumed to be radial symmetric, 
it can be written as B z (r,B rnax ) = B rnax f(r). The 
true field is then simply B = B z (r, B max ) cos{(3). 
The error is then given by 



A new vacuum chamber which allows to cover a 
angular spread of ±45° or —10° to 80° is being de- 
signed. 



3.2 Measurement methods 

To check the field homogeneity the magnetic field of 
the spectrometer was mapped for different coil cur- 
rents and different positions. To measure the field's 
componentes a transversal and a longitudinal hall 
probe (± lmT) were used. 



AB =( cos 2^(^ A r) 2 + 
+cos 2 f3 (gf^ AS) 2 + B 2 Z sin 2 (3 (A/3) 2 )s 

Using gfgt- = f(r) and g£g « (3 2 the relative 
error is than given by 



For the experiment this becomes 
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Therefore the uncertainty depends 
strongly on the position in the field. At the mag- 
net edge the relative error is up to ~ 15%. For 
the magnet's bulk it is on the other hand in the range 
1°/ 
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Since the magnet's poles are quite parallel and 
have a large diameter (~ 170 mm) compared to the 
pol gap (~ 25 mm) it likely that the field in the mag- 
net's bulk does not vary much in z and has only small 
radial component. Therefore it is more important to 
map the field at the magnet's edge. For simplicity the 
field has only been measured in the (r,z) plane. The 
stray field measurements have been done after remov- 
ing the original vacuum chamber. To determine the 
exact position of the hall probe the probe was fixed 
to a table which was able to be adjusted in height 
with ±lmm accuracy. The radial position was con- 
trolled by a measuring stick which was glued to the 
table (±lmm). The experimental setup is given in 
Fig. [4j The radial field was measured by using a lon- 
gitudinal and the vertical field by using a transversal 
hall probe. 



Since the radial field can be very small and the 
vertical field is large (« 1.150 T) a small off zero an- 
gle between the z axis and the hall probes surface 
gives a fkrx that can not be neglected. Therefore the 
(O^e^urjpd radial field has to be corrected by this flux 
contributed from the vertical field. The contribution 
from the vertical field was determined by inserting 
the probe inside the magnet's bulk where the radial 
field should be zero. This allows to calculate a simple 
correction term. 

3.3 Results on Magnetic properties 

The magnetic field perpendicular to the poles B z 
has been measured as a function of the radial po- 
sition Fig [5] If the magnetic field is normalized by 

f = B "B^tT h ' {t can be seen ( Fi g® that the field's 
normalized shape does not change with the coil cur- 
rent. 




Figure 4: Experimental setup to map the B z ans B r 
component of the magnetic stray field as a function 
of r and z. The uncertainty in r and z is 1 mm. 




Figure 5: The magnetic field as distance from the 
magnet's centre. The measurements were done along 
the connection line between port 4 and 8 




Figure 6: The results from keeps its shape even for 
low coil current. 



This allows to introduce the magnetic radius as 



magnetic 



f(x) dx 



with B max as the maximal field for a certain coil 
current. Due to the measurement the magnetic radius 
Rmacnetic = (167.0 ± 0.6) mm can be found. 

The length of the stray field AR str ay defined as 
the region where the field drops from 95% to 5% was 
determined as AR str ay = (110 ± 5) mm. The ra- 
tio between the stray field length and the magnetic 
radius was found to be ^ Rstra « = 0.66 ± 0.02 

To make sure that the field is radial symmetric the 
measurement from Fig. [5] has been done for different 
directions. The radial symmetry can clearly be seen in 
Fig. [3 
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Figure 7: The magnetic field as a function of the dis- 
tance to the magnet's centre for different directions 
as showed in the inserted sketch. The measurement 
have been done for 47 A coil current. 



Far from the edge inside the magnet the field is 
parallel to the z axis and far outside the field is negli- 
gible. Therefore the field has only been measured at 
the magnet's edge as a function of its radial and ver- 
tical position. The results can be seen in Fig. [5J The 
magnetic field lines are quite parallel. The measured 
radial field is smaller than 15% of the vertical field. 
As it will be seen in the part on simulation correc- 
tions due to the vertical position and the radial field 
are small. 




Figure 8: The absolute value of the magnetic field 
as a function of the radial and the vertical position. 
The vertical position z — is the midpoint of the pol 
gap. The measurement have been done at 87 A coil 
current. 

The magnetic field for the spectrometer was 
simulated in Poisson/Supcrfish [5 . The pro- 
gram Poisson/Superfish solves the two dimensional 
Maxwell equation for certain boundary conditions. 
The boundary conditions were in such a way chosen 
that they represent a (x,y) cut though the magnet's 
centre. The resulting field lines are given in [9] Since 
the coil consists of 140 windings a total current of 
12180 A = 140 • 87 A was used for the calculations. 
A maximal field of 1.1 5T could be found. This is in 
good agrement with the (1.150 ± 1)T, which where 
found in the measurements. 




Figure 9: Simulated field for the two dimensional 
Maxwell equations, with the boundary conditions 
given by a x-y cut through the zero point. 

To determine the ramping characteristic and check 
on memory effects the magnetic field was changed adi- 
abaticly. To make sure that changes were adiabatic 



the coil current was slowly (5Amin~' 1 ) shifted. Two 
cycle starting at OA, going to 47 A, 87A respectively 
and back to OA were measured. The magnetic field 
for the larger cycle is given in Fig. [TO] It can be found 
that 

B{I 11 {t en d)) < B(I l2 (t enc i)) 

with / 7l (t) < I l2 {t)\/t < t end and I 7l (t end ) = 
I~f 2 (tend)- 7 is the path that was used to change the 
current. The differences B{I l2 {t end )) — B{I 11 {t end )) 
are shown in Fig. [TTJ 



4 Discussion on Simulation 

4.1 Simulation methods 

To determine the trajectory of an electron in the spec- 
trometer it is required to take the real field of a mag- 
net into account. It was assumed that the magnetic 
field is constant for each iteration step. Using the 
relativistic momentum p = mo 7 v the equation of 
motion is given by 

dp 

dt= qVAB 
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Figure 10: The plot shows the magnetic field in the 
magnets centre vs. the coil current (0 Ato87 A). It 
can bee seen that the field is nearly proportional to 
the current. 



Figure 11: The plot shows the difference for the 
B{I l2 (t en d)) — B(I 11 (t en d)) as a function of the coil 
current. The current was changed slowly (5 Amin~ ) 
to stay in the adiabatic limit. 
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Since there is only a Lorentz force acting on the 
electron 7 and (3 remain constant. Using this the dif- 
ferential equation simplifies to 



™ 7 (7 2 /3 2 + l) 0- 
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cPx 
ds 2 



= qf3cn v A B 



Where the chain rule has been applied to write 
the problem as a function of the trajectory length s. 

With this assumption the equation of motion can 
be solved for each iteration step. To get an idea of the 
algorithm's robustness an electron's trajectory was 
simulated in a homogeneous magnetic field. Since an 
analytic expression for an electron's motion in such 
a field exists (p = ^j) it is possible to compare the 
derivation between the analytical and numerical so- 
lution Fig. Q2] It can be seen that for reasonable 
step size the true value of p is achieved with a high 
accuracy (w 1% ). 




Figure 12: To estimate the simulation's accuracy the 
trajectory of an electron in a homogeneous field has 
been calculated and compared with the theoretical 
value. Where s is the step size, p is the theoretical 
bending radius and p S im is the simulated. 



For all the simulation an experimental setup as 
shown in Fig[l3] was used. In the simulation either 
the magnetic field measured by Danfysik or the one 
measured by us were used. 
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Figure 13: Schematically picture of the setup used for 
the simulation. This setup have been chosen for the 
all the simulations. The detector is in a cylindrical 
shape placed behind the spectrometer. The magnet's 
physical diameter is 540mm the magnetic diameter is 
330mm. 



Figure 14: The trajectories for electron beams at dif- 
ferent energies in the real magnetic field (B ~ 1.15 T). 
The broken circle indicates the spectrometer's mag- 
netic radius. 
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Figure 15: The plot shows the energy E vs. the 
bending angle a. 



4.2 Results from simulation 



For the measured magnetic field the deflection of an 
electron beam can be simulated as a function of its 
energy as seen in Fig. Q3] If this trajectories are 
compared with the trajectories for a field given by 



B — Bn 



flection angle of maximally 2% can be found Fig 



an error in the de- 
101 



Since an error of 2% is to large the problem can not 
be handled analytically (as in the appendix) but has 
to be treated numerically. All the further discussions 
will base on numerical simulation. 



Figure 16: The simulated results (cross) and the an- 
alytical calculation (line) are in a good contradiction 



(±2%). 

This enables us to plot the bending angle vs. the 
electron's energy (Fig. 1 1 5(1 - For large bending radius 
p 3> Rmagnetic the analytic expression as derived in 
the appendix can be simplified. 

a = arccos y p p2+R2 J = arccos I - | ^ ' y j 
Taylor 

This shows that at large energies the bending an- 
gle a is inverse proportional to the energy. Therefore 
the resolution decreases dramatically with increasing 
energy. 
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Figure 17: The plot shows the energy resolution 
vs. the electron beam's energy E. A focusing effect 
can be seen at around 110 MeV. For a setup with 1 
meter drift length between source and magnet edge 
and an additional drift length of 0.2 meter between 
the magnet edge and the screen. The magnetic field 
was the real measured field at 87 A (1.15 T). 



4.2.1 Simple energy resolution 

A mono energetic beam with a certain divergence 
spread passes a drift length and then the spectrom- 
eter. During the propagation the beam diverges and 
the signal in the detection plane gets spread. For 
a free drift and if space charge effects are neglected 
the spread after a distance Z is given by A = Z A'. 
Therefore an increase in divergence results in a de- 
crease in energy resolution. This can be seen by a 
simple case of a none mono energetic beam. If a beam 
with a spread in divergence passes the experimental 
setup, the beam will be transformed as indicated in 
Fig. [2TJ While the beam passes the dipole magnet 
electrons with higher momentum get deflected less 
than electron with smaller momentum ( p = j . 
Therefore the electrons are separated in space accord- 
ing to there momentum. This is indicated in Fig. 1211 
This allows to calculate the spectrometer's reso- 
lution in a simple way. To each point y in the de- 
tection plane a certain energy range AE can be as- 
signed (green line in Fig. l2Tj) . The energy resolution 
is than defined as . For a setup with 1 meter drift 
length between source and magnet edge and a addi- 
tional drift length of 0.2 meter between magnet edge 
and detector the energy resolution vs. the energy can 
be calculates for different divergence (as given in Fig. 
I17j) . Since the spread in y is proportional to the drift 
length the energy resolution depends crucial on the 
distance between source and magnet (Fig. 




Figure 18: Sketch for three different path. A is a 
electron going through the magnets centre, B and C 
cut the magnetic area off centre due to there finite 
divergence. It can be seen that B is staying in the 
magnetic region for the shortest time, where C is in 
the magnetic region for the longest time. The differ- 
ence in path results in a focusing. 

It can be seen in Fig. [T7|that at a certain energy 
a focusing effect occurs. This can be illustrated by 
looking on two rays, one of the ray is following the 
central line A and the other ray is drifting to the left 
site B, due to its finite divergence. Since the magnet 
has circular shaped pole the pathway in the magnet 
for the ray A is longer than the pathway for B (see 
Fig. [18]) . Therefore ray A gets more deflected than B. 
On the other hand ray C which is drifting to the right 
of the central ray is staying in the magnetic field for 
a longer time than A and B and gets therefore more 
deflected. This is similar to the effect described in [3] . 



A E/E 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

D . h . „ ' I m l 

maget-deteclor 



Figure 19: Resolution vs. beam energy and distance 
between spectrometer and screen. The setup uses a 
real magnetic field generated by 87 A coil current and 
1 meter drift length between source and magnet edge. 
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Figure 20: Resolution vs. beam energy and distance 
between spectrometer and screen. The setup uses a 
real magnetic field generated by 87 A coil current and 
1 meter drift length between source and magnet edge. 



It can be seen in Fig. [19] how the focusing effect 
is depending on the drift length between magnet and 
screen. A focusing is achieved mainly on energies be- 
tween 100 MeV and 200 MeV. The results shown in 
Fig. [19] have been made without taking space charge 
effects into account [3]. 
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Figure 21: A schematic picture of the evolution of a 
area in p-y space after passing the spectrometer. 

4.2.2 Lifting the degeneracy 

However if the divergence is known the degener- 
acy can be lifted. Measuring the divergence can 
be achieved by putting a screen Ti in front of the 
spectrometer (Fig. |2"2"1) . This allows to measure the 
beam's size and position. Now the opportunity is 
given to simulate the propagation of different mono 
energetic electron beam's (at energy {Ei}i e j) with 
the determined divergence and initial position. Once 
the propagation of mono energetic electron beam is 
known the distribution S(Ei) on the detection screen 
Ti can be calculated easily. Since it is known how 
mono energetic beams map on the detection screen Ti 
as a function of their energies the opportunity is given 
to write an arbitrary signal S, which is detected on 
T 2 , as a linear combination of mono energetic beams 
S(Ei). 

S = A l S(Ei) 

S can be fitted to the measured distribution S by 
a root mean square method, with the restriction that 
Xi > for all i € I. The set {Xi}i^i is than equal the 
beams energy distribution. 




Figure 22: The sketch gives a proposal for an ex- 
perimental setup to measure the energy resolution 
more accurate. T\ and Ti illustrate two screens which 
measure the beam's position and spatial distribution 



simultaneously. For simplicity the picture shows a 
beam consisting of electrons at two different energies. 
After a certain drift length D\ the beam has due to 
its divergence a certain spatial distribution. After 
passing the spectrometer the electrons get deflected 
in a known way according to their energy. The result 
measured on the screen is a superposition of electrons 
from different energies. 

To get a rough understanding what influence noise 
at the detection screen T2 on this analyse procedure 
has, a mono energetic beam E with a uniform distri- 
bution in divergence (from —3 mrad to 3 rarad) was 
studied. The noise was simulated by broadening the 
signal in y direction by 0.5 mm. The width of the 
electron energy distribution was determined and uses 
for calculating the energy resolution (Fig. 
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Figure 23: 

This results are only very rough estimation. The 
inaccuracy from the screen T\ has not been taken in to 
account. Also only mono energetic beams were stud- 
ied. However this maybe far from real experimental 
data. 



4.2.3 Motion in z direction 

The bulk field does not have any influence on the mo- 
tion's z component. Therefore only the stray field 
accelerates the particle along the z direction. Due to 
[3] the stray field at the magnet edge causes a defo- 
cusing in the z direction. Indeed this can be observed 
in the numerical simulation [2H The defocusing as a 
function of the energy can be simulated (Fig. 




Figure 24: The x component vs. the z compo- 
nent of the trajectory. It can be seen that the elec- 
tron propagates freely except at the magnet edge 
where it gets bended. For the simulation no space 
charge effect have been included. The simulated elec- 
tron had none divergence and a initial position at 
(1.27, —0.001, 0)) to. For the spectrometer the real 
magnetic field at 87 A coil current was used. 
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Figure 25: The displacement in z direction yield- 
ing from bending effects at the magnet's edge vs. 
the beam energy. For the simulation no space 
charge effect have been included. The simulated elec- 
tron had none divergence and a initial position at 
(1.27, —0.001, 0)) to. For the spectrometer the real 
magnetic field at 87 A coil current was used. 



However the broadening created by the magnet 
edge is small compared with the beam's divergence 
that can be expected. 



5 Conclusion 

For the magnets physical properties it can be stated 
that the magnetic radius is (167. 0± 0.6) mm. A stray 
field occurs over a length of (110 ± 5) mm. However 
the magnetic field's shape f(r) does not depend on 
the coil current. Furthermore the maximal magnetic 
field B z was measured for different coil current and 
it was shown that the field is proportional to the cur- 
rent. In addition a memory effect for the magnetic 
field could be observed, depending on the path that 
was chosen to adjust a certain current. In the adi- 
abatic limes the different between lowering and ris- 
ing direction was found to be smaller than 0.015 T. 
For further studies it could be of interest to measure 
memory effects beyond the adiabatic limit. 

Numerical simulation have shown how the bend- 
ing angle depends on the electrons energy. It can 
be concluded that the bending angle does not change 
much for large energies (E > 1 GeV) . This decreases 
the energy resolution dramatically for increasing en- 
ergies. However at low energies the resolution is rel- 
atively high (l°/ 00 for E = 150 MeV). Since some 
interests lies in electron energies up to 5 GeV a pro- 
posal has been made to resolve energy spread in this 
regime energies. However much more studies have to 
be done to clarify if this procedure is appropriated or 
not. 

In addition focusing effects in y direction and defo- 
cusing effects in z direction has been studied. It was 
shown that this effects are stronger on low energies 
(< 300 MeV). 

6 Appendix 

If it is assumed that the magnetic field is homoge- 
neous in a circular area with radius Rmagnetic and 
zero other wise (B = B max e z X^ x 2 +y 2 <RmagnM ), a 
analytical solution for the bending angle can be found. 
For simplicity it is further assumed that the incident 
beam is perpendicular to the magnet's surface. Sym- 
metry reasons request that the outgoing beam is per- 
pendicular as well. Since the magnetic field in the cir- 
cular region is constant the electron's trajectory will 
be a segment of a circle. The bending radius can eas- 
ily be calculated by p = ^L. Knowing the bending 
radius, the magnetic radius and the angle of incident 
{(3 = 0) allows us to calculate the bending angle a 
(see FiglH- 



P 2 -R 2 



The simulated results for the bending angle are 
plotted on a scale which shows a deviation from the 
analytical value as a deviation from the straight line 
Fig. HH Even though that the analytical method 
gives quite good results the accuracy is not high 
enough. 




Figure 26: The sketch shows the trajectory of 
an electron in a magnetic field given by B = 
B e z x r~2~i — 5> o ■ Inside the magnetic field the 

electron follows a circle and outside the electron fol- 
lows a straight line. 
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